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Technological advances in recent years in the control of 
quantum coherence have supposed an important stimulus 
to the field of 'quantum chaos', namely, the study of the 
effect of classical chaos in quantum mechanics. A corner- 
stone result in quantum chaos is the so called Bohigas- 
Giannoni-Schmit (BGS) conjecture [H that states quan- 
tum spectra of classically chaotic systems are universally 
correlated according to the Wigner-Dyson random ma- 
trix ensembles (WD) Q for scales comparable to the 
mean level spacing, A. Similarly, according to the Berry- 
Tabor-Gutzwiller conjecture [H , the level statistics of sys- 
tems whose classical counterpart is integrable are well 
described by Poisson statistics, namely, the spectrum is 
uncorrelated. 

Despite its unquestionable success, these conjectures 
are notjdways verified. Exceptions include the Harper 
model 0], and the ubiquitous kicked rotor In the 
latter, for short time scales, both quantum and classical 
motion are diffusive in momentum space. However, quan- 
tum diffusion is eventually arrested due to destructive in- 
terference effects that localize eigenstates in momentum 
space. In this limit spectral correlations are described 
by Poisson not by WD statistics. This counter-intuitive 
feature is usually referred to as dynamical localization 

SI. 

Deviations from the BGS conjecture are also expected 
for eigenvalue separations 5E ~ h/ts due to weak local- 
ization effects The typical scale tE = A -1 1 log h\ is the 
Ehrenfest time with A the classical Lyapunov exponent. 
However, these corrections, though universal, do not re- 
ally invalidate the BGS conjecture since 8E ~ h/tE ^> A. 

From the previous discussion it is clear that the main 
deviation from the semiclassical picture contained in the 
BGS conjecture is due to localization effects. This fact 
rises the following questions: Why in some quantum 



chaotic systems localization is only weak but in others 
is strong enough to induce a metal-insulator transition? 
Is it possible to give a more precise relation between clas- 
sical motion and quantum features? Are there any other 
universality classes in the context of non-random Hamil- 
tonians? The aim of this paper is to address these ques- 
tions within the framework of the one parameter scaling 
theory [§] introduced originally in the context of disor- 
dered systems [§]. 

A key concept in this theory is the dimensionless con- 
ductance g introduced by Thouless [icj . It is defined 
either as i) the sensitivity of a given quantum spectrum 
to a change of boundary conditions in units of the mean 
level spacing A oc 1/L d , or ii) as g = E c /A where, E c , 
the Thouless energy, is an energy scale related to the 
diffusion time to cross the sample. In the metallic limit 
E c = fiD c i as /L 2 (D c i as is the classical diffusion constant) 
and therefore g oc L d ~ 2 . On the other hand if the particle 
is exponentially localized due to destructive interference, 
g oc e~ L /£ where £ is the localization length and L is the 
system linear size. 

The change of g(L) with the system size, (3(g) — 
d ^dfoL i P rov ides information about localization. For 
L — > oc, (3(g) = d — 2>0ina metal (without quantum 
corrections) and (3(g) = log(<?) < in an insulator. The 
OPT is based on the following two simple assumptions: 
i) The (3(g) function is continuous and monotonous; ii) 
The change in the conductance with the system size only 
depends on the conductance itself. 

With this simple input the OPT is able to predict 
three universality classes in the L — > oo limit: For d > 2 
and disorder sufficiently weak, g — > oo. The system has 
metallic like features: Eigcnfunctions arc dclocalized in 
real space, and the spectrum is correlated according to 
WD statistics. For d < 2, or d > 2 and strong disorder, 
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g — > 0. The system is an insulator: Eigenfunctions are 
exponentially localized and the spectrum is correlated ac- 
cording to Poisson statistics. For d > 2, a metal-insulator 
transition takes place at a certain density of impurities 
and energy. It is characterized by a size independent 
dimensionless conductance g — g c , namely, (3(g c ) = 0. 
Eigenstates at the transition are multifractals and the 
spectral correlations are universal [I~0 ] , but different from 
WD and Poisson statistics. We now study how the OPT 
can be adapted to the context of quantum chaos. 

One parameter scaling theory in quantum chaos. -The 
first problem is the very definition of the Thouless energy. 
In disordered systems it is estimated after an average 
over many disorder realizations. In quantum chaos no 
such ensemble average can be carried out. However, a 
Thouless energy can still be defined provided: i) The 
ensemble average is replaced by an average over initial 
conditions; ii) The classical phase space is homogeneous 
with no islands of stability. This condition guarantees 
that the Thouless energy does not depend on the initial 
conditions chosen. 

The localization problem in quantum chaos is defined 
in momentum space, not in real space. This can be un- 
derstood as follows: In a classically integrable systems 
the number of conserved quantities (canonical momenta) 
is equal to the dimensionality of the system. In quantum 
mechanics each of these canonical momenta becomes a 
good quantum number which labels the state. An in- 
tegrable system is thus localized in momentum space in 
the sense that there exists a basis of momentum eigen- 
states in which the Hamiltonian is diagonal and conse- 
quently the spectrum is uncorrelated (Poisson statistics). 
As classical symmetries are reduced, the Hamiltonian is 
no longer diagonal in any momentum basis and even- 
tually a transition from Poisson to WD statistics is ex- 
pected. 

After these clarifications we are ready to define g in 
non-random Hamiltonians. As the classical dynamics is 
by no means restricted to standard diffusion, we study 
the more general case (p 2 ) ~ i M ((...) stands for aver- 
age over initial conditions) with p, > 0. The Thouless 
energy is given by E c oc N^ 2 ^ 1 where N is the system 
size, namely, the size of the basis of momenta eigenstates 
in which we express the Hamiltonian. The mean level 
spacing is in many cases given by A ~ 1/N d , but there 
are important exceptions: i) for periodic potentials, the 
Bloch theorem applies, the spectrum is continuous, and 
A = 0; ii) for systems whose eigenstates are exponen- 
tially localized, A / even in the N — > oo limit; iii) 
in systems with a singular continuous spectrum the scal- 
ing with the system size may be anomalous A oc N~ d / d ' 
(d e ^ 1). A precise definition of d e may depend on the 
system in question. In the Harper model [4], d e rs 1/2 
stands for the Hausdorff dimension of the spectrum. Wc 
are now ready to define the dimensionless conductance in 
quantum chaos. In cases i) and ii) above, g — > oo (metal) 



and g = (insulator) respectively. In case iii) (including 
d r = 1), 
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The running of g is thus described by (3(g) — ^p^jf = 
Idas — j- — Under the assumptions of the OPT 
and using Eq. ([T]) we propose the following alternative 
definition of universality class in quantum chaos: i) If 
Idas > 0, eigenfunctions are delocalized as in a metal and 
the spectral correlations are described by WD statistics, 
ii) If "/das < 0, eigenfunctions are localized as in a insula- 
tor and the spectral correlations are described by Poisson 
statistics, iii) If ^das — 0, eigenfunctions are multifractal, 
a metal insulator transition takes places and the spectral 
correlations are universally described by critical statistics 
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Several remarks are in order: i) Universality is re- 
stricted to scales of the order of the mean level spacing; 

ii) If d e 1, or for a continuous spectrum A = 0, it 
may not be possible to carry out a statistical analysis 
of the spectral correlations. In these cases our classifi- 
catory scheme still holds but the characterization of the 
system as a metal or an insulator must rely directly on 
the eigenfunctions statistics or the transport properties; 

iii) Quantum destructive interference may modify g. For 
instance, all eigenstates are localized in a 2d disordered 
system despite the fact that semiclassically g oc L d ~ 2 is 
constant as at the Anderson transition; iv) In order to 
define the Thouless energy uniquely the moments of the 
classical distribution of probability must be described by 
a single scale. 

Stability of the semiclassical predictions.- The previ- 
ous discussion highlights the need of a clear understand- 
ing of the effect of Anderson localization on (3(g). Wc 
focus on the d — 1 and d e = 1 case, the general- 
ization d > 1 is straightforward. We investigate the 
quantum dynamics of a classical diffusion process de- 
scribed by the following fractional Fokker-Planck equa- 
tion 0, - D clas £^) P(p,t) = S(p)S(t) where 

D c i as is the classical diffusion coefficient and p is a real 
number. The moments of the distribution P(p, t) (if 
well defined) are {\p\ k ) oc t fJ " k / 2 , with k a real num- 
ber. The classical propagator in Fourier space (t <->■ w, 
P q) is given by K (q,uj = 0) = Dci J g 2/ M where 
v stands for the spectral density. One-loop corrections 
due to interference effects [l3[ take the form K~ 1 (q) = 



K^(q)~ 



|fc| 2 /M 



The quantum diffu- 



l o w 2 

sion coefficient D quant to this order is easily obtained by 
performing the integral above: 
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FIG. 1: (Color online) (a) Comparison of the classical and 
quantum (|p| fe ) oc t^=( fc ' for the case of a Id kicked rotor 
with V(q) oc log \q\, \q\ a ■ In agreement with the analytical re- 
sults the classical time dependence is not modified by quan- 
tum corrections for a < 0. (b) Spectral rigidity A3 (L) for 
V(q) = 10 log |g| . The level statistics in this case have all 
the signatures of a metal insulator transition such as scale 
invariance and linear As(L) ~ with x ~ 0.33 < 1. 



where C is a different constant for each case. The im- 
portance of the quantum effects depends strongly on the 
value of 2/ 11. 

In the region 1 < 2/fi < 2, quantum corrections di- 
minish the value of the classical diffusion constant. A 
renormalization group analysis shows that the semiclas- 
sical prediction 0(g) < for L — > 00 still holds but the 
transition to localization is faster due to quantum cor- 
rections. In the region < 2/fi < 1 quantum corrections 
are subleading with respect to the classical term. As a 
consequence the semiclassical prediction /3(g) > is not 
altered by quantum localization effects. For 2//i = 1, 
the logarithmic behavior resembles superficially that of 
2d disordered system, D quant = D c i as — Cln(L/Z) where 
exponential localization occurs. However there are dif- 
ferences. For anomalous diffusion the mean free path, I 
in the 2d case, is the inverse momentum q" 1 . Therefore, 
the corrections to the bare coupling constant are small for 
small momentum q ~ 1/L. Qualitatively this implies the 
absence of eigenstate localization. For a rigorous proof 
based on the evaluation of (3(g) including higher order 
terms we refer to 131 and references therein. Therefore 



the quantum properties of an Id Hamiltonian with clas- 
sical dynamics such that (p k ) oc t k resemble those of a 
disordered conductor at the metal-insulator transition. 
In summary, the semiclassical predictions of a Id system 
whose classical dynamic is anomalous are stable under 
quantum corrections. In higher dimensions this might 
change. For instance, in d > 2, we expect deviations due 
to localization for j c i as > 0. 

Examples. - We now test the predictions of the OPT 
in different deterministic systems. 

Id Kicked rotor with classical singularities. - We study 

2 

the quantum dynamics of, W = \ + V(q) S(t — nT) 
for different non-analytic potentials lj,|l5|, V(q) — e\q\ a 
and V(q) — elog(|g|) with q S [—it, it), a e [—1,1] and e 
a real number. The classical evolution is dictated by the 



map p n+1 
InRef. [I 



Pn 



dq„ 



q n +i = q n + Tp n+ i (mod 2ir). 
it was found that for a > 0.5 classical diffu- 
sion is normal. Therefore (p 2 ) oc t and jcias = — 1 < 0. 
Thus the OPT predicts Poisson statistics even though 
the classical dynamics is chaotic. It is remarkable that 
by using scaling arguments dynamical localization can be 
predicted without having to map the problem onto a Id 
Anderson model. 

For —0.5 < a < 0.5, classical diffusion is anomalous, 
(pk) ^ t Ki- a )_ Using Eq [JJ 



1 -a 
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For a > the OPT predicts Poisson statistics no matter 
what the effect of quantum corrections is. For a < we 
expect WD statistics since quantum corrections do not 
modify qualitatively (3(g). For a = (log singularity), 
we expect an Anderson transition. Therefore dynamical 
localization can be overcome, even in one dimension, if 
classical diffusion is fast enough; i.e. \i > 2. 

These predictions have been tested by studying the 
quantum evolution operator U over a period T in a basis 



of plane waves \n), (m\U\n) 



J_ e -i2-rrMn 2 /N 
N 



where (f>(l,m,n) = 2n(l + 8 )(m - n)/N - iV(2n(l + 
9 )/N), I = -(N - l)/2, . . . , (N - l)/2 and < 6 < 1; 
8q is a parameter depending on the boundary conditions 
(9q = for periodic boundary conditions). The eigenval- 
ues and eigenvectors of U are computed by using standard 
diagonalization techniques. In Fig. la it is observed that, 
in agreement with the perturbative analysis, the classi- 
cal time dependence of (p k ) is not modified by quantum 
corrections for a < 0. Thus a genuine Anderson transi- 
tion is expected for a = 0. The scale invariance of the 
spectrum and the analysis of the level statistics (see Fig. 
lb) fully confirms the theoretical prediction. 

3d Kicked rotor.- We study the quantum dynam- 
ics of a 3d kicked rotor with a smooth poten- 
tial (the 2d version was studied in Ref. [13] ) : W = 
\(txp\ + T 2 p\ + T 3 p 2 3 ) + V(qi,q 2 , q 3 ) E„ <*(* " nT ) with 
V(qi,q2,qz) = k cos(qi) cos(q 2 ) cos(q 3 ) and Ti,T 2 ,t 2 in- 
commensurate. The spectrum of the evolution matrix 
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FIG. 2: (Color online) (a) Quantum (p 2 ) in the 3d kicked ro- 
tor with a smooth potential for three different kicking values: 
above, below and at the transition k = k c ~ 3.3. In agree- 
ment with the OPT (p 2 ) oc t 2/3 at k = k c . (b) The spectral 
rigidity As(L) versus L for the same kicking values. A transi- 
tion from WD to Poisson statistics is clearly observed as the 
kicking strength decreases. The level statistics at k ~ k c have 
all the signatures of a 3d Anderson transition such as scale 
invariance and Aa(L) ~ xL/15 with x ~ 0.19 < 1. 



was obtained by evolving a quantum state IV'(O)) an d 
performing the Fourier transform of (ip(t)\ip(0)). Classi- 
cally the diffusion (p 2 ) oc t is normal provided that the 
classical phase space is fully chaotic. Quantum dynamics 
depends strongly on A:. In analogy with a 3<i disordered 
system, we expect destructive interference stop the clas- 
sical diffusion for sufficiently small k. In the opposite 
limit, quantum effects are small and diffusion persists. A 
careful finite size scaling analysis 12] has confirmed this 
picture [161 ]. We have found a metal-insulator transition 
at k = k c « 3.3. According to the OPT, since (3(g) = 
at the transition, quantum diffusion must be anomalous 
(p 2 ) oc t 2 / 3 . Likewise, level statistics are described by 
WD (Poisson) statistics in the limits k ^> (<C)A; C . For 
k = k c it is expected spectral correlations be similar to 
those of a 3d disordered system at the transition. As is 
shown in Fig. 2, the numerical results fully agree with 
these theoretical predictions. 

Harper model.- The Id Harper model 0], H = cos(p) + 
A cos(27rer2;), with a irrational, undergoes a metal in- 
sulator transition at A = 2. Classically the system is 
integrable. However, the quantum motion is diffusive 



(a; 2 ) ~ t 2dH with dn the Hausdorff spectral dimension. 
With this information and Eq. (QJ we can compute the 
dimensionless conductance and (3(g). As was expected 
(3(g) = and g — g c is size independent as it is expected 
at the metal insulator transition. Our simple method 
predicts correctly the metal insulator transition in this 
model as well. In conclusion we have investigated under 
what conditions the one parameter scaling theory can be 
utilized in quantum chaos. We have utilized it to de- 
termine the number of universality classes in quantum 
chaos and propose a more accurate definition of them. 
The universality class related to the metal insulator tran- 
sition has been investigated in detail. We have tested 
our theoretical predictions in different kicked rotors and 
the Harper model. Our findings open the possibility of 
studying the metal insulator transition experimentally in 
a much broader class of systems. 

AMG acknowledges financial support from a Marie 
Curie Outgoing Action, contract MOIF-CT-2005-007300. 
AMG thanks Igor Rozhkov for illuminating discussions 
and a critical reading of the manuscript. JW is grateful to 
Prof. C.-H. Lai for his encouragement and support, and 
acknowledges support from Defence Science and Tech- 
nology Agency (DSTA) of Singapore under agreement of 
POD0613356. 



[1] 

m 

[3] 

[4] 
[5] 

[6] 

[7] 
[8] 

[9] 
[10] 

[11] 
[12] 

[13] 

[14] 

[15] 

[16] 
[17] 



0. Bohigas, M.J. Gianonni, C. Schmit, Phys. Rev. Lett. 
52, 1 (1984). 

M.L. Mehta, Random Matrices, 2nd Ed., Academic Press 
(San Diego, 1991). 

Chaos in Classical and Quantum Mechanics, M. C. 
Gutzwiller, Springer- Verlag, New York (1990); M.V. 
Berry, M. Tabor, Proc. R. Soc. London A 356, 375 
(1977). 

M. Ya. Azbel, Zh. Eksp. Teor. Fiz. 46, 929 (1964); D. R. 
Hofstadter, Phys. Rev. B 14, 2239 (1976). 
G. Casati, B. V. Chirikov et al, in Stochastic Behavior in 
Classical and Quantum Hamiltonian Systems, edited by 
G. Casati and J. Ford, Lecture Noted in Physics Vol. 93 
(Springer, Berlin, 1979). 

S. Fishman, D. R. Grempel, R. E. Prange, Phys. Rev. 
Lett. 49, 509 (1982). 

1. L. Aleiner, A. I. Larkin, Phys. Rev. B 54, 14423 (1996). 
E. Abrahams, P.W. Anderson et al, Phys. Rev. Lett. 42, 
673 (1979). 

P.W. Anderson, Phys. Rev. 109, 1492 (1958). 

J.T. Edwards, D.J. Thouless, J. Phys. C 5, 807 (1972). 

V.E. Kravtsov, K.A. Muttalib, Phys. Rev. Lett. 79, 1913 

(1997); S. Nishigaki, Phys. Rev. E 59, 2853 (1999). 

B.I. Shklovskii, B. Shapiro, et.al. Phys. Rev. B 47, 11487 

(1993). 

A. D. Mirlin, Phys. Rep. 326, 259 (2000); A. D. Mirlin, 
Y. Fyodorov et al, Phys. Rev. E 54, 3221 (1996). 
J. Liu, W. T. Cheng, C.G. Cheng, Comm. Theor. Phys. 
33, 15 (2000). 

A.M. Garcia- Garcia, J. Wang, Phys. Rev. Lett. 94, 
244102 (2005). 

A.M. Garci'a-Garci'a, J. Wang, In preparation. 

E. Doron, S. Fishman, Phys. Rev. Lett. 60, 867 (1988). 



